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We consider the embedding of the supersymmetric Standard Model with broken R-parity in the 
minimal supergravity (mSUGRA) model. We restrict ourselves to the case of broken lepton number, 
the B3 mSUGRA model. We first study in detail how the tree-level neutrino mass depends on the 
mSUGRA parameters. We find in particular a strong dependence on the trilinear supersymmetry 
breaking A-parameter, even in the vicinity of the mSUGRA SPSla point. We then reinvestigate 
the bounds on the trilinear R-parity violating couplings at the unification scale from the low-energy 
neutrino masses including dominant one-loop contributions. These bounds were previously shown to 
be very strict, as low as 0(10^'') for SPSla. We show that these bounds are significantly weakened 
when considering the full mSUGRA parameter space. In particular the ratio between the tree-level 
and 1-loop neutrino masses is reduced such that it may agree with the observed neutrino mass 
hierarchy. We discuss in detail how and in which parameter regions this effect arises. 
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I. INTRODUCTION 

The experimental observation of neutrino oscilla- 
tions, and thus of neutrino masses, is an experimental 
indication that the Standard Model of particle physics 
(SM) is incomplete [lM3- 

Experimentally, neutrinos must be relatively light. 
Direct laboratory measurements restrict their masses 
to be below 0(10 MeV - leV) 0411, depending on 
the flavor. Cosmological observations even give up- 
per bounds of CfOTeV) on the sum of the neu- 
trino masses [3, [l2| - [l4j . Furthermore, the atmospheric 
and solar neutrino oscillation data are best fit if the 
squared neutrino mass differences are C'(10~^eV^) and 
O(10~^eV^), respectively 0, [H]- This allows for one 
massless neutrino. 

In principle, it is easy to extend the SM Lagrangian 
by a Dirac neutrino mass term [7|. However, right- 
handed neutrinos and new Yukawa couplings of 0{^ 
10~^^) are in this case needed. Such tiny couplings 
seem to be very unnatural and might point towards 
a dynamical mechanism, that explains the small neu- 
trino masses. Furthermore, the right-handed neutri- 
nos can have an unspecified Majorana neutrino mass. 

Most prominently discussed are extensions of the 
SM involving the see-saw mechanism, by introducing 
right-handed neutrinos and fixing the new Majorana 
neutrino mass scale to be large, cf. Refs. jTl [l6l - l2l| . 
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The see-saw mechanism is also naturally incorporated 
into supersymmetry (SUSY) mm. 

Supersymmetry is one of the most promising ex- 
tensions of the SM. It is the unique extension of 
the Lorentz spacetime syrnmetry, when allowing for 
graded Lie algebras [23, 1231 ■ Furthermore, it provides 
a solution to the hierarchy problem of the SM J26| - |3C| . 
More importantly here: neutrino masses can be gen- 
erated without introducing right-handed neutrinos if 
lepton number is violated, cf. for example Refs. |3ll - 

sii. 

The most general gauge invariant and renormaliz- 
able superpotential of the supersymmetric extension 
of the SM with minimal particle content (SS M) p os- 
sesses lepton number conserving (LNC) terms [42, \^ 
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and also lepton number violating (LNV) terms 
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where i,j,k = 1,2,3 are generation indices. We have 
employed the standard notation of Ref. [4J| . 

The LNV interactions violate the discrete sym- 
metries R-parity and proton-hexality (Pg), however, 
they conserve baryon triality (B3) |45l - l48l |. Note that 
B3 stabilizes the proton because it suppresses the 
baryon number violating interactions. R-parity, Pe 
and B3 are the only discrete symmetries, which can be 
written as a remnant of a broken anomaly free gauge 
symmetry |45| - |48| . In the following, we assume that 



B3 is conserved and thus R-parity and Pg are violated. 
Eq. (IlJ and Eq. ([2]) constitute the full renormalizable 
superpotential allowed by this symmetry. For reviews 
of such theories see for example Refs. |49l - (5l| . 

Beside the superpotential, also the soft-breaking La- 
grangian of the B3 conser ving SSM exhibits lepton 
number violating operators [52| 
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where again i, j, fc — 1,2,3 arc generation indices. L, 
£, Q and D are the scalar components of the lepton 
doublet, lepton singlet, quark doublet and down quark 
singlet superfield, respectively. Furthermore, hu (hd) 
denotes the up-type (down-type) scalar Higgs field. 
Beside the term proportional to mj = , the operators 

in Eq. (jSj arc the soft-breaking analog of the terms in 
Eq. ©. 

The LNV terms in Eq. ^ and Eq. dS]) lead to the dy- 
namical generation of neutrino masses. For example, 
the bilinear terms in Eq. ^ mix the Higgsinos, the su- 
persymmetric partners of the Higgs bosons, with the 
neutrino fields and thus generate one non-vanishing 
neutrino mass at tree-level |3ll - [33l |37H40| . 

In this paper, we derive bounds on the trilinear LNV 
couplings of the superpotential, Eq. ^, from the up- 
per cosmological bound on the sum of neutrino masses 



A^- 



<0.40eV, 



(4) 



at 99.9% confidence level. The bound was deter- 
mined by a combination of the Wilkinson Microwave 
Anisotropy Probe (WMAP) and Large Scale Structure 
(LSS) data. 

In order to perform a systematic study, we restrict 
ourselves to the well motivated framework of the B3 
minimal supergravity model (mSUGRA) |52| . which 
provides simple boundary conditions for the SSM pa- 
rameters at the grand unification scale (Mqut)- We 
describe the model in the next section in detail. We 
employ the full set of renormalization group equations 
(RGEs) at one loop |52H5a | in order to obtain the B3 
SSM spectrum and the neutrino masses at the clcc- 
troweak scale (Mew)- We then derive bounds on the 
LNV trilinear couplings at Mqut- 

Bounds on trilinear LNV couplings within this 
model were also derived in Ref. [52| from the genera- 
tion of neutrino masses at tree-level. It was claimed 
that neutrino masses put an upper bound of 0(10"^ — 
10~^) on most of the trilinear couplings in Eq. ([2]). 



However, it was shown in Ref. [5j| that the tree-level 
neutrino mass can vanish in certain regions of the B3 
mSUGRA parameter space. In our analysis, we espe- 
cially focus on these regions of parameter space. We 



show that the bounds on the trilinear couplings can be 
weakened up to O(10~^), depending on the boundary 
conditions. 

Wc go beyond the former work in several aspects. 
Beside the tree-level neutrino mass, we also include 
the dominant contributions to the neutrino mass ma- 
trix at one-loop. These contributions were neither in- 
cluded in the calculation of the bounds in Ref. [52[ nor 
in Ref. [SJ]. However, as we show in Sec. llVi the loops 
dominate in the regions of parameter space where the 
tree-level mass vanishes. They must thus be included 
when determining the bounds. 

In Ref. [5JI there is only a brief explanation of the 
dominant effect that leads to a vanishing tree-level 
mass in B3 mSUGRA. We give for the first time a de- 
tailed and complete explanation of how different con- 
figurations of the B3 mSUGRA parameters at Mqut 
can affect the tree-level and loop contributions to the 
neutrino masses at A/ew ■ Although we restrict oursclf 
to the framework of B3 mSUGRA, the mechanisms de- 
scribed in this publication also work in more general 
models. Furthermore, we calculate bounds for all tri- 
linear LNV couplings, whereas Ref. |5J| focused only 



on the couplings Ai33 and A'gg. We also update the 
bounds given in Ref. [53 according to the more recent 
and stronger bound on the sum of neutrino masses, cf. 
Eq. dl. 

Going beyond the work presented here, we believe 
our results can help find LNV SUSY scenarios that 
explain the observed neutrino masses and mixing an- 
gles. Within the framework of B3 mSUGRA, Ref. [HI 
searched for a minimal set of LNV parameters which 
can explain the measured neutrino parameters. They 
found sets of five parameters [two trilinear LNV cou- 
plings together with the three mixing angles that de- 
scribe the lepton Yukawa matrix, cf. Eq. Ip] that give 
the right masses and mixing angles. Ref. [41| claimed 
that the tree-level mass is always much larger than the 
loop induced masses. But we show in the following, 
that the loops can exceed the tree-level masses in B3 
mSUGRA. Therefore, it should be possible to find a 
smaller set of LNV parameters that lies in this region 
of parameter space and thus posses much larger LNV 
couplings than those found in Ref. |41|. However, an 
investigation of the complete neutrino sector is beyond 
the scope of this paper and will be postponed to a fu- 
ture publication. 

We finally note that (large) trilinear LNV couplings 
can lead to distinct collider signatures at the Large 
Hadron Colhdcr (LHC), e.g. 

• Supersymmetric particles (sparticles) can be 
produced singly at a collider, possibly on res- 
onance j5ll |56|463| . For example, single resonant 
slepton production at the LHC via A^^^, Eq. ^ 
di-Ill 113, HI. An excess over the SM back- 
grounds is visible if X'--f. > 0(10"'^), depending 
also on the sparticle masses [53, [S^, ISj, IfiJ ■ 

• A LNV coupling A^fc (A^^) of > ©(IQ-^) at 



Mqjjt can significantly change the running of 
the sparticle masses, such that the scalar elec- 
tron or muon (sneutrino) is the LSP [52|, [65l - l68j . 
This can dramatically change the SUSY collider 
signatures, because (heavy) sparticles normally 
cascade decay down to the LSP [6^ SM ■ 

Note that the generation of neutrino masses via the 
bilinear terms in Eq. ([2]) and the corresponding collider 
signatures have also been investigated; see for example 
Refs. [70-|81| and references therein. 

This paper is organized as follows. In Sec. |lll we 
review the parts of the B3 mSUGRA model that are 
relevant for this work. Sec. IIIII then shows the differ- 
ent contributions to the neutrino mass matrix that we 
employ to derive the bounds. We explain the main 
mechanism leading to a vanishing tree-level neutrino 
mass in B3 mSUGRA in Scc. lIVI and derive the bounds 
on the LNV trilinear couplings in Sec. |Vl These sec- 
tions are the central part of our paper. We conclude 
in Sec. |Vll 

App. |X] explains the additional subleading depen- 
dence of the neutrino masses on the B3 mSUGRA pa- 
rameter not described in Sec. HVl 



II. THE B3 MSUGRA MODEL 

The general B3 SSM has more than 200 free param- 
eters [S^l ■ This large number is intractable for detailed 
phenomenological studies. For that purpose the sim- 
plifying B3 mSUGRA model was proposed in Ref. [5^ , 
which we now discuss. 



A. Free Parameters 

In the B3 mSUGRA model the boundary conditions 
at Mqut arc described by the six parameters 

Mo, Ml/2, Ao, tan/3, sgn(^), A, (5) 

with 

Ae{A,,fe,A^J. (6) 



Here Mg, M^^ and Aq arc the universal scalar mass, 
the universal gaugino mass and the universal trilinear 
scalar coupling at the grand unification scale (A/gut), 
respectively, tan /? denotes the ratio of the Higgs vac- 
uum expectation values (vevs) w„ and Vd, and sgn(/i) 
fixes the sign of the bilinear Higgs mixing parameter 
fi. The magnitude of fi is determined dynamically by 
radiative electroweak symmetry breaking (REWSB) 
[83| . These five parameters are the conventional free 
parameters of the R-parity or proton- hexality conserv- 
ing mSUGRA model H. 

In order to incorporate the effects of the LNV inter- 
actions in Eq. ^ and Eq. ([3]) exactly one additional 
non-vanishing trilinear coupling A e {Aijfe,A^j,} is 



assumed at A/qut- Further LNV couplings are gen- 
erated via the RGEs at the lower scale. Note, that 
the bilinear couplings Ki and Di arc both set to zero 
at Mqut via a basis transformation of the lepton and 
Higgs superfields [3l[ . (For the most general case of a 
complex rotation see Ref. [83 ■) This is natural for uni- 
versal SUSY breaking [S^]- However, at lower scales ki 
and Di are generated via the RGEs [33]; see Sec. IIIGI 
The complete low energy spectrum is obtained 
by running the RGEs down from Mgut to 
A/ew- For that purpose we employ the program 
SOFTSUSY-3.0.12 j.86, 87]. We calculate the neutrino 
masses with our own program. Note that we work in 
the CP-conserving limit throughout this paper. 



B. Benchmark Scenarios for Parameter Scans 

We center our analysis around the following B3 
mSUGRA parameter points 

Point I: A/1/2 = 500 GcV, A/q = 100 GeV, 



tan^ = 20, sgn(/i) = +1, ^0 = 900 GeV, A = A: 



233 



Point II: Ml/2 = 500 GeV, Mq = 100 GeV, 
tan^ = 20, sgn(/i) = +1, Aq = 200 GeV, A = A233 

Point II differs from Point I only by the choice of 
the LNV coupling and the size of Aq. We have chosen 
these points as examples because the tree-level con- 
tribution to the neutrino mass is small around Point 
I and II and therefore one-loop contributions are im- 
portant. Both points lead to squark masses of 0(1 
TeV) and slepton masses of around 300 GeV, with a 
scalar tau (stau) as the LSP. 

Note that in the LNV SSM a stau LSP is as well mo- 
tivated as a neutralino LSP [52, 58, 65|, IMS] ■ Either 
will decay via the LNV interactions and cosmological 
constraints do not apply j91| . 

In addition, we ensured that both points lie in re- 
gions of parameter space where various other exper- 
imental constraints are fulfilled, such as the lower 
bound on the lightest Higgs mass from LEP2 [92, IS^ 
and constraints from the anomalous magnetic mo- 
ment of the muon [9J], from 6 — >■ 57 [95|, and from 
Bs -^ A*^M^ [Sal; see Sec.|V]for details. 



C. Renormalization Group Equations and 
Radiative Electroweak Symmetry Breaking 

An important feature of the B3 mSUGRA model is 
that lepton number violation leads to mixing between 
the lepton superfields Li and the Higgs superfield Hd- 
Furthermore, sneutrinos, the superpartners of the neu- 
trinos, can acquire vevs Vi (i = 1,2,3). Note that it 
is possible to rotate away the Ki terms in the super- 
potential at any given energy scale by an orthogonal 
rotation of the fields £„ = {Hd, Li) [3l|, [52|, l85i] . 



The corresponding bilinear soft-breaking terms pro- 
portional to Di, Eq. ([3]), can be rotated away in con- 
junction with Ki if Di and Ki are aligned. This con- 
dition is fulfilled at Mqut in the B3 mSUGRA model 
if the underlying supergravity superpotential satisfies 
the quite natural condition [52] 



f{Z'i;ya) = fl{Zt) + fliVa) , 



(7) 



where the superfields Zi belong to the observable sector 
and the superfields ya to the hidden sector. 

However, when evolving the parameters down to the 
weak scale, Ki, D^ ^ are generated via the RGEs. 
The leading terms for A € {A^j.} are given by [53 | 



dt 
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and 
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Here t = ln((5//io) with Q the rcnormalization scale 
and /io an arbitrary reference scale, /i^ j, = Aq x A^ ^. 

at Mqut, c/. Eq. ((3)). _B is the soft supersymme- 
try breaking analog of the Higgs mixing parameter ^ 
and (h(7)33 [(h£;)33] is the soft-breaking analog of the 
Yukawa coupling {Yu)33 [(Yb)33] J52|. gi and 52 {Mi 
and M2) are the U(l)y and SU(2) gauge couplings 
(gaugino masses), respectively. We see in Eqs. ^ and 
^ that the RGEs differ, and therefore Ki and Di will 
no longer be aligned at the weak scale [33|. The case 
As {Ayfc} is analogous up to the color factor 3. 

The sneutrino vevs Vi, the bilinear Higgs parameter 
|/z| and the corresponding soft breaking term B are 
determined by REWSB, which has been discussed in 
detail in Ref. |52| for the LNV case. 



Neglecting higher order corrections |96l - l98| . which 
are n ot im portant for the following qualitative discus- 
sion lllol. the sneutrino vevs can be written as 15211 



where {rn.f^'^) is the squared soft-breaking lepton dou- 
blet mass matrix and q = \/3/5 Oi . m? - originates 

from the LNV soft-breaking Lagrangian, Eq. ([3]). It 
mixes the down-type Higgs fields, hj^, with the lepton 
doublet scalars, Li, and is zero at Mqut- That is, be- 
cause we take within mSUGRA the mass matrix of the 
fields Ca = (hd, Li) to be diagonal and proportional 



to Mq at Mgut- However, m^ - 7^ is subsequently 
generated via the RGEs, cf. Eq. ([55)1 . 

As we will see in Sec. IHI 11 sneutrino vevs and non- 
zero bilinears Ki lead to neutrino masses at tree-level 
because they mix neutrinos and neutralinos. 



D. Quark Mixing 

The RGE evolution of the parameters in the B3 
mSUGRA model from Mgut to A/ew depends on the 
Higgs- Yukawa coupling matrices Y^;, Y^ and Yjj, 
cf. Eqs. dll) and ^. In particular, the RGEs of the 
LNV violating parameters are coupled via the non- 
diagonal matrix elements of the Higgs- Yukawa cou- 
plings. Therefore a knowledge of the latter is crucial 
for the analysis of bounds on the LNV parameters. 

The initial parameter set of the B3 mSUGRA model 
at A/qut is given in the electroweak basis so that for 
the RGE evolution the Higgs- Yukawa couplings (or 
the quark- and Icpton-mass matrices) are also needed 
in the electroweak basis. However, from experiment 
we only know the masses and the CKM matrix 

VcKM = Ui Dl (12) 

at Mew- Here Ul (Dl) rotate the left-handed up- 
(down-) quark fields from the mass eigenstatc basis to 
the electroweak basis. For simplicity, we take Y^i and 
Y[/ to be real and symmetric and thus the rotation 
matrices for the right-handed quark fields are identical 
to the ones for left-handed quark fields, Ur ~ Ul 
and Dr = Dl . Because of the uncertainty about the 
neutrino masses and mixings we will assume a diagonal 
Ye in the weak basis. 

When determining the neutrino masses, we will con- 
sider two limiting cases at M^w, following Ref. jSalSSl 
HI: 

• "up— type mixing" the quark mixing is only in 
the up-quark sector, 

Ul,R = VcATA/, Dl.r = 1 , 



(Mihv, 



m 



hdLi 



yUKj 



Vd + DiVu , (10) 



with 



{M?h 



(mL^)y + K,Kj + -Ml cos 2(3 5,^ 



(5'+g|)„. 2 



sin' 



^pY.'^u.,, (11) 



Yu X Ud = 
Y 



Vd = d\a.g{md,ms,mb) , 



jj xvu^ VcKM ■ diag{mu,mc, rrit) ■ Vckm ■ 



(13) 



• "dovirn— type mixing" the mixing is only in the 
down-quark sector. 



D 



L,R 



V, 



CKM, 



u 



L.R 



Y_D x Vd = ^CKM ■ diag(md, ms,mb) ■ Vckm 1 



Yu y. Vu = diag(m„, nic, rrit) 



(14) 



Here md,ms,m{, (rnu,mc,mt) denote the masses of 
the down-type (up~type) quarks. 

The choice between up- and down-type mixing has 
a strong effect on the final results for the LNV cou- 
plings A S {^ijk} with j ^ k, as we will show in SeclVl 
(see Tab. HI . The reason is that the generated tree level 
neutrino mass is proportional to the off-diagonal ma- 
trix element (Y^))^^,, cf. the discussion in Scc. lIIII and 
Sec. lIVI Our results (for the tree-level neutrino mass) 
in Sec. |V] can be easily translated to scenarios which 
lie between the limiting cases of Eqs. ([T5|) and (ITil) . 
One only needs to know the respective Yukawa matrix 
elements {YD)jk- 



III. NEUTRINO MASSES 





In this paper, we investigate bounds on lepton- 
number violating couplings at Mqut within the B3 
mSUGRA model, which arise from the generation of 
too large neutrino masses at Mew- We therefore need 
to identify the dominant contributions to the neutrino 
masses. 

It was stated in Ref. [S^] that the main contribu- 
tion stems from mixing between neutralinos and neu- 
trinos, which leads to one non- vanishing neutrino mass 
at tree-level, cf. Sec. IIII 11 However, as we will show 
in the next two sections, this is only true in parts of 
the B3 mSUGRA parameter space. It is possible that 
the different terms in the tree-level mass formula can- 
cel each other. We then need to identify the dominant 
contributions, which arise at one-loop. 

A complete list of all one-loop contributions is given 
in Ref. [38i] . where they are formulated in a basis- 
independent manner. Most of the one-loop contribu- 
tions are proportional to the mass insertions that mix 
the neutrinos with the neutralinos. They thus also 
vanish when the tree-level neutrino mass vanishes and 
are negligible in the region we are interested in. 

The remaining dominant one-loop contributions are 
on the one hand due to loops involving two R-parity 
violating vertices and are thus either proportional to 
A^ or to A'2, cf. Fi g. [J fTTTJ . We will review these 
contributions in Sec. IIII 21 On the other hand, loops 
with virtual neutral scalars (i.e. Higgscs and sncutri- 
nos) and neutralinos, which are shown in Fig. [2l can 
also give large contributions to neutrino masses. These 
loops are proportional to the mass difference between 
CP-even and CP-odd sneutrinos, cf. Sec. IIII 31 

According to Ref. [S^l, there is in principle also a 
contribution which is proportional to A x 1?^. However, 
this contribution is suppressed by two ore more orders 
of magnitude in the regions of parameter space where 
the loops dominate over the tree-level mass. Note that 
Di vanishes near the tree-level mass minimum as we 
will show in Sec. IIVI We therefore neglect it in the 
following. 

Further one-loop contributions arc only present in a 
lepton- and Higgs-supcrfield basis with non- vanishing 



FIG. 1: Loop contributions to the neutrino mass matrix 
via a non-vanishing product of B3 couplings Xikn x \jnk 
(upper figure) and A^j.^ x Aj„fc (lower figure). See Sec. IIII 2] 
for more details. 



HI 




Al 




xl 



FIG. 2: Loop contributions to the neutrino mass matrix via 
a non-exact cancellation of loops with GP-even and GP- 
odd neutral scalars. Note, that there is a relative minus 
sign between the two diagrams. See Sec. IIII 31 for more 
details. 



sneutrino vevs. This is the case in the B3 mSUGRA 
model. However, we have checked that in our param- 
eter scans, these contributions are at least one order 
of magnitude smaller than the dominant one and thus 
negligible for calculating the bounds. Note, that they 
are also aligned with the tree-level mass, because the 
sneutrino vevs vanish near the tree-level mass mini- 
mum, cf. Sec. IIVAI 

We conclude, that the contributions to the neutrino 
masses which we review in the following, are sufhcient 
to calculate the correct bounds on the LNV couplings 
A and A'. However, in order to calculate the correct 



neutrino mass spectrum and mixing angles, all one- 
loop contributions given in Ref. [38| must be taken 
into account. This lies beyond the scope of this paper. 



1. Tree-Level Contributions 

In the context of the B3 mSUGRA Model, neutrino 
masses are generated at tree-level through mixing be- 



tween neutrinos and neutralinos. Analogously to the 
standard see-saw mechanism 0, [l6l - l2ll | (with the neu- 
tralinos taking over the role of the right-handed neu- 
trinos) , an effective 3x3 neutrino mass matrix is gen- 
erated [nisi. 



J 



M-:„ = 



^l{^hgl + A'hg^) 



2vuVd{Migi + M252) - 2^MiM2 



AiAi A1A2 A1A3 
A2A1 A2A2 A2A3 
A3A1 A3A2 A3A3 



(15) 



where Mi (M2) is the bino (wino) soft-breaking mass 
and 



Ai = Vi -Vd — , 



i = l,2,3. 



(16) 



This matrix has one non-zero eigenvalue which can at 
Mew be simplified to j53 | 



167raGUT J2i=i ^f 



Ml 



/2 



(17) 



if we take into account the gaugino universality as- 
sumption at A/gut, leading to M2 = I^Mi = 

-^Mi/2 at Mew IH. Here acuT = 5GUT/47r « 
0.041 is the grand unified gauge coupling constant [52l |. 



2. Contributions from \X- and X A -Loops 

In the region of parameter space where the tree- 
level neutrino mass, Eq. (|17p . vanishes, loop induced 
neutrino masses give the dominant contributions. As 
we will show in Sec. |Vl the most important loops are 
those proportional to the product of two LNV trilin- 
ear couplings. The corresponding squark-quark and 
slepton-lcpton loops arc shown in F ig. FT] The result- 
ing neutrino mass contributions are |35| 



im',% 



327r2^ 

k.n 



Kkn^onkmi^ sin2(^^ln 



327r- 



■ X! Kkn>^']nkmdk sin 24)i In 



k^n 




where m^^ i'nidk ) f-re the lepton (down-quark) masses 
of generation k, and 0^ i'Pn) the mixing angles that de- 
scribe the rotation of the left- and right-handed slep- 
ton (down-squark) current eigenstates of generation n 



to the two mass eigenstates, m^ and m^ {nij^ and 
TO J ), respectively. Note that the squared sfermion 
masses are linear functions of the mSU GRA parame- 
ters Mq and M^,2J see for example Ref. jlOOf . For the 
calculation of Eq. ((T8| and all following calculations, 
we have used the two- component spinor formalism as 
described in Ref. [loH . 

For the first two sfermion generations, the sfermion 
mixing angles are small and we approximate Eq. (|18p 
by using the mass insertion approximation (MIA) as 
described in Ref. [41|. The slepton (and down-squark) 
mass eigenstates are replaced by the respective Icft- 
and right-handed eigenstates with mass mi and 
771J . The mixing angle can be approximated by 



sin 2(f>'^ 



2{Mf 



where 



LR.\2 



{Mn 



^ir. 



(hfi). 



/itan/3 



(19) 



(20) 



denotes the left-right mixing matrix element of the 
charged sleptons of generation n. {h.E)nn is the trilin- 
ear soft-breaking analog of the lepton Yukawa matrix 
element (Ye)„„ [H. 

A similar formula is obtained for sin 205^. One only 
needs to replace in Eq. ((TO)) and Eq. (^0)) I -^ d, I -^ 
d, (Yb)„„ o (Y_d)„„ and (h£;)„„ o (h_D)„„, where 
(h£))„„ is the soft-breaking analog of the down-quark 
Yukawa matrix element (Yi3)„„. 



3. Contributions from Neutral Scalar-Neutralino-Loops 

As the final source of neutrino masses, we consider 
contributions arising from loops with neutral scalars 
and neutralinos, cf. Rcfs. jSJ-jSg. Most important is 



the contribution from sneutrino-antisneutrino mixing, 
as we will see in Eq. (|26l). 

If CP is conserved, sneutrinos i>i and antisneutrinos 
v* mix to form CP-invariant mass eigenstates 



.--,+ 



1 






(21) 
(22) 



,--.± 



If lepton number is conserved, the i'j masses arc 
degenerate and the CP-even (CPE) and CP-odd 
(CPO) contributions to the neutrino mass from neu- 
tral scalar-neutralino-loops cancel, cf. Fig. [5] 

In contrast, if lepton number is violated, the v,^ 
masses are in general different, so the cancellation is no 
longer exact. This is due to the fact that the CPE and 
CPO neutrinos mix differently with the CPE and CPO 
Higgs fields, respectively. The size of this contribution 
to the neutrino masses is roughly proportional to the 
mass splitting Atti?. = 771?+ — 77i?_, cf. Eq. ([^5]) and 

Refs. Hm. 

The neutral scalar-neutralino-loops, shown in 
Fig. [21 lead to the following contributions to the neu- 
trino mass matrix kill 



{<% 



1 ^ ^ 



327r' 



k=l L=l 



-^(2+j)L^(2+j)lBo(0, m^o , to|o ) 

'(23) 

where m^o (fc = 1 . . . 4) are the neutralino masses and 
A^ is the 4x4 neutralino mixing matrix in the bino, 
wino, Higgsino basis |l02| . The two-point Pass arino - 
Veltman function is conventionally denoted Bq |l03l |. 
m^o (to^o ) with L = 1, . . . , 5 are the mass eigenval- 
ues of the CPE (CPO) neutral Higgs bosons and CPE 
(CPO) sneutrino fields. They can be obtained with 
the help of the unitary matrix Z~^ {Z~), which diago- 
nalizcs the mass matrices of the CPE (CPO) neutral 
scalars, i.e. 



{Z+fM 



CPE 



Z+ = 



= dia.g{mlo,mjjo,ml+,ml+,ml+) 



diag(m^o ) 



(24) 



and 



(z- 



fMcpoZ- 



diag(m^o , m\o , m^- , to^- , m^_ ) 



diag(TO^o ) 



(25) 



see Ref. [4l| for additional details. 

In order to analyze the dependence of this contribu- 
tion on the mSUGRA parameters, we make use of the 



fact that in the BamSUGRA model, Eq. ^ can be 
approximated by [39[ 



(<%, 






fe=i 
Am? 



■In 



(™1 -TO^o)' 



by expanding around 7717,0 and 777 



^L>2 



(26) 



The mass 



splitting, Ar77?., in Eq. ([26]) between CPE and CPO 
sneutrinos of generation i is then given by [36| 



A777? 



-AB^M^ml sin^ /3 



(777^0 " mljimlo ~ 777? J (777^0 " m?J 
^^ (Bv, - b.Vdf 



{vl + vl){B^+D^) 



(27) 



4. Numerical Implementation 

The numerical calculation of the neutrino mass ma- 
trix was done in the following way. We first employed 
SOFTSUSY-3 . . 12 pj , H^j to obtain the low energy 
mass spectrum |118l | . We then used our own program 
to calculate the neutrino mass matrix. The tree-level 
contribution was derived from Eq. (|15p . For the AA- 
and A'A'-loops, we employed Eq. p^ . if third gener- 
ation sfermions were involved. However, for sfermions 
of the first two generations we used the MIA as given 
in Eqs. (fTO l) and ([20]) . 

For the neutral scalar-neutralino-loops, we in prin- 
ciple employed Eq. ((23)) . However, instead of per- 
forming the large numerical cancellation between CPE 
and CPO neutral scalars directly [square bracket in 
Eq. (1231)] I we used an MIA to calculate the devia- 
tion from exact cancellation in the R-parity conserving 
(RPC) limit, following Ref. [41|. The resulting formula 
is quite lengthy and we refer the interested reader to 
Ref. [41| for details. We have cross checked our pro- 
gram with the help of Eq. ((261) and Eq. ^^. All our 
calculations are performed in the CP-conserving limit. 



IV. !^-MASSES: DEPENDENCE ON MSUGRA 
PARAMETERS 



In the literature it has frequently been assumed 
that the tree-level contribution to the neutrino mass, 
Eq. (fTT)) . in the B3 mSUGRA model dominates over 
the loop contributions, cf. for example Refs. |4ll.l52|. 
However, as has been noted in Ref. [5J], in certain re- 
gions of B3 mSUGRA parameter space, the tree-level 
neutrino mass vanishes even when k^ ^ 0. 

We demonstrate this effect in Fig. [31 where we dis- 
play the tree-level neutrino mass (solid red line) as 
a function of Aq. The other B3 mSUGRA parame- 
ters are given by Point I with A233IGUT = 10~^, cf. 



Sec. Ill BI We see that the tree-level mass, m^J'^'^, van- 
ishes around ^o ~ 910 GeV. In the vicinity of this 
minimum, m^J°'^ drops by several orders of magnitude 
over a wide range of ^Oi ^nd it is therefore not a (large) 
fine-tuning effect. In this case the loop contributions 
will dominate the neutrino mass matrix, resulting in 
much weaker bounds on the involved A coupling, cf. 
Sec. |Vl Thus the bound crucially depends on the 
choice of Aq. 

We emphasize that the range of Aq for which weaker 
bounds may be obtained is quite large. In an interval 
of AAq « 100 GeV around the minimum, we obtain 
bounds on A233 that are at least one order of magni- 
tude smaller than the bound derived at for exainple 
Aq = GeV. Much weaker bounds can therefore be 
obtained without a lot of fine tuning. 

In this section, we aim to explain in detail the origin 
of this cancellation, considering as an explicit exam- 
ple the case A g {A^-^.}. We focus on the dependence 
of ■m*IJ'^° on the mSUGRA parameter Aq, because it is 
always possible to find a value of ^0 [for a given set of 
parameters tan/3, M1/2, Mq, and sgn(ju)] such that the 
tree-level neutrino mass vanishes. All arguments can 
analogously be applied to a Xijk coupling, as discussed 
in App. IA5I Note for the further discussion that we 
can always obtain a positive A by absorbing a possible 
sign of A via a re-definition L ^ —L and E ^ —E oi 
the lepton doublet and lepton singlet superfields, re- 
spectively. We also note that the generated neutrino 
masses scale roughly with A^, cf. the following dis- 
cussion. 



A. ^0 Dependence of the Tree— Level Neutrino 
Mass 
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FIG. 3: Ao dependence of m^J"'' and the terms determining 
the sneutrino vev V2, Eq. (|31|l . at the REWSB scale (used 
in SDFTSUSY to calculate the sneutrino vev). Since the 
scale affects the parameters only logarithmically, there are 
only minor changes when running to A/ew- The other B3 
mSUGRA parameters are that of Point 1 with A233IGUT = 
10"^ SecHTH 



\vi\ , i — 2 on Aq. It possesses a clear minimum which 
is close to the minimum of to|;''°°. 

This behavior can be understood by taking a look 
at the (tree-level) formula for the vev Vi, Eq. pO)). For 
A G {'^iyfe} it can be written as 



Wc now discuss the dependence of the tree-level 
neutrino mass at Mew as a function of Aq at Mqut- 
Recall from Sec. lIII ij that 



1 



m„ 



ocA^ 



Vt - Vd- 



M 



(28) 



From the RGE of Ki, Eq. (|S]), we obtain as the domi- 
nant contribution 

K, a /iA^fe {Yn),k ^ ^A^, ^^^^ (29) 

Vd 

at all energy scales, where {md)jk denotes a matrix 
element of the down quark mass matrix. Therefore, 



■yd— ex A-j-fe • {md)jk , 
A* 



(30) 



without further dependence on mSUGRA parameters. 
Thus, the dependence of the tree-level neutrino 
mass, Eq. (^5]) . on the mSUGR A pa rameters is solely 
through the sneutrino vev Vi |ll9| . In Fig. [31 the 
dashed green line explicitly shows the dependence of 



iM?h 



biVu - (m^ =: -I- ^lKi)Vd 



with 



(M?),, = (niL^),, + -M| cos 2/3. 



(31) 



(32) 



Here, we have neglected terms proportional to k| and 
vf , because they are much smaller than (m.^'^)ii and 
Mf . Note that we only obtain one non-zero sneutrino 
vev because A'- -^ violates only one lepton flavor. 

In many regions of parameter space the sneutrino 
vev in Eq. (pS)) is at least two orders of magnitude 
larger than the term w^Kj/zi. Thus the minimum of the 
neutrino mass can only occur when the sneutrino vev 
is drastically reduced. As we shall see, the sneutrino 
vev becomes very small, when there is a cancellation 
between the two terms in Eq. ([5T|) . 

The second term of Wi in Eq. ([5T|) . (m^ - + ^Ki)vd, 

and the prcfactor \/{Mf)ii are always positive and 
depend only weakly on Aq. This can be seen in Fig. |3] 
for (m^ - +fiKi)vd (dotted-dashed blue line) and also 



for l/{Mf)ii (solid turquoise line). This behavior can 

be easily understood: 

The soft breaking parameter, mj = , Eq. (j3]), is zero 

I— -I 
at Mgut and is generated at lower scales via [52| 






dt 



(33) 

where J-{fh?) is a linear function of the soft-breaking 
scalar masses squared and of the down-type Higgs 
mass parameter squared, /i^-^. [(h£))j7j] is the soft- 
breaking analog of A'^-^. [(Yd)^^] with /i^^-j, = X^^j, x Aq 
[(h.D)jk = O^D)jk X Aq] at A/gut- The second term 
in Eq. ([33]) thus depends on Ag. However, J^{fh?') is 
in general much larger than A^ due to several con- 
tributions from soft breaking masses [521. Therefore, 
varying Aq does not significantly change the magni- 
tude of m\ - as long as ^o is not much larger than 
the sfermion masses. 

Concerning the term ^ki in (mj - -I- finAvd, we 
note from the RGE for Kj, Eq. ([S]), that the only Aq 
dependence of k^ stems from its pro port ionality to fi. 
fi at Mew can be approximated by [100l | 
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FIG. 4: Running of /1233 for various values of Aq. The other 
B3 mSUGRA parameters are that of Point I, Sec. Ill Bl with 
A233IGUT = 10"^ and Ml/2 = 500 GeV. 



^^ 



ciM^ 



-C2Mi%- 



C3Al + CiAoMi/2-^ . (34) 



Here ci and C2 are number s of C (l) whereas C3 and C4 
are only of 0(10^1-10^2) [H^. Therefore, except for 
Aq ^ A/q, A/1/2, the order of magnitude of fi remains 
constant when varying Aq. 

We conclude that (mj - -f- finAvd depends only 

weakly on Aq and therefore, Di is decisive for the 
Aq dependence of the vev Vi and thus of 7nJ;''°°. If 
the first term in Eq. pip . DiVu, is positive and only 
slightly larger than the (nearly constant) second term, 
(m^ - -t- ^Ki)vd, Vi can equal VdK.i/^ and we get 

to';'=^ = 0, cf. Eq. dill). 

The strong ^0 dependence of the magnitude of 
DiVu is also displayed in Fig. [3] (dotted magenta 
line). We observe that lUiWul is often larger than 
(m^ - + fiKi)vd (dotted-dashed blue line). However, 
near the tree-level neutrino mass minimum (solid red 
line), it drops below (mr - + fiKi)vd and Vi can equal 
Vd — - In this case m^J°'^, Eq. ((28|) . vanishes. 

In order to understand this behavior of Di, we need 
to understand how Di is generated via the RGEs. Re- 
call that Di = a,t Mgvt within the B3 mSUGRA 
model. The generation of Di primarily depends on 
the running of the trilincar soft breaking mass /i^ j. 



,dD 



16^^'-^ - ^6ti{YD),kK 



dt 



ijk 



(35) 



We find the contribution i n Eg . © proportional to B 
is typically much smaller |12l| and we here focus on 



the effects due to h'. 



The dominant terms of the 
corresponding RGE are given by [52|, |66| 



ijk- 



dt 



16 



gi {2M3K.k ~ h') + 



(36) 



where g^ {M3) denotes the SU(3) gauge coupling 
(gaugino mass). At A/qut this equation simplifies to 



167r 



2dKjk 

dt 



y5GUT(2Mi/2-Ao)A^fe + .... (37) 



Keeping for now all paramete rs ex cept Aq fixed 



(with sgn(^) = +1 and A-^j. > |122| ). we can clas- 
sify the running of h'--f^ , Eq. ([55)) and Eq. ([57)) , in the 
following way (see also Ref. [g^ for a detailed discus- 
sion): 



(a) Ao ^ 2M1/2 (including negative values oi Aq): 
Since the right hand side (RHS) of the RGE for 
h'^ji^, Eq. ([36)1 . is always positive and large, h'^j/, 
is quickly reduced from its initial value of Aq x 
Kjk ^^"-^ even becomes negative when running 
to lower energies. This behavior is displayed in 
Fig- ID (dashed green line), where the running of 
/1233 is shown for different boundary conditions 
at AfouT- 



(b) Ao ~ 2M1/2: If the size of ^0 is comparable 
to 2Mi/2, h'^jf. will be fairly constant at high 
energies, cf. the dotted magenta line in Fig. |3| 
However, when running to lower energies it will 
still start decreasing, but more slowly than in 
case (a). This is due to the fact that Afa and 
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FIG. 5: Running of the bilinear coupling D2, Eq. ((35}, for 
the same parameter sets as those in Fig. U 



Kjk themselves increase significantly (by factors 
of approx. 2.5 and 3, respectively; see Ref. [6g) 
when running to lower energies. Thus the term 
2M3X'^j^. eventually dominates in Eq. ([36]) even 
if initially ^0 ^ 2Ai^i/2- This leads to a small, 
negative h'^-j, at low energies. 

(c) Ao > 2M1/2: h[-^ is large at A'/qut and is fur- 
ther increased when running to lower energies. 
This is due to the negative RHS of the RGE for 
^iife' Elq. ([55)) : see also the dotted-dashcd blue 
line in Fig. S) 

Caveat: Since the term 2M^\[jf. in Eq. ([36)) in- 
creases by a factor of approximately 8 « 3 • 2.5 
when running from A/qut to Mew [as men- 
tioned in (b)], h[-i, only strictly displays the be- 
havior of case (c) when Aq > 2OA/1/2. Other- 



wise, h'--^. will decrease once the term 2My,\'^-f. 
dominates. 

Because Di is zero at Mqut and, according to 
Eq. ([55)) . also proportional to the integral of h[-f, over 
ln((5), points (a) 
for A: 



(c) have the following consequences 



(a) Ao <C 2M1/2: Since h'^-j, always becomes neg- 
ative below some energy scale close to A/qut, 
the RHS of Eq. ()35p is positive. This leads to a 
large negative Di at Mz as can be seen in Fig.O 
(dashed green line). Consequently, all terms ex- 
cept DiVu become negligible in Vi, Eq. pip , and 
thus \vi\ at Mew is large, dominating the tree- 
level neutrino mass, Eq. ([25)) . 

(b) Ao w 2M1/2: Due to the initially negative RHS 



of Eq. ([35)) at energies close to A'/qut (where 

^'iik ~ ^0 X •^i,fe)' ^i ^'^s* increases when run- 
ning to lower energies but then starts decreas- 
ing once h[-f, becomes negative, c/. the dotted 
magenta lines in Fig. |3] and Fig. |5| At some 
energy scale Q, Di becomes small such that Vi, 
Eq. ([?!)) . can equal w^ — . A cancellation between 

these two terms in to|^''°°, Eq. ([25)) . at the scale Q 
will then occur. This corresponds to a vanishing 
tree-level neutrino mass if Q ~ Mew- 

(c) Aq > 2M1/2: The RHS of Eq. ([35|) is always 
negative with a large magnitude such that we 
get a large positive Di at the weak scale, c/. the 
dotted-dashed blue line in Fig. |5j As in case 
(a), DiVu provides the main contribution to \vi\, 
Eq. ([5T)) . Therefore, \vi\ is large and dominates 



TO, 



Eq. ^. 



Summarizing, the tree-level neutrino mass has a 
minimum in the parameter region where the size of 
Aq is comparable to 2Mi/2. This is mainly due to the 
running of the parameters Di and h'^^^. that affect the 
sneutrino vevs; in particular due to a partial cancel- 
lation in Eq. ([55)) . Note that in Fig. [3] the tree-level 
neutrino mass vanishes at Aq w 910 GeV, which is 
indeed close to 2Mi/2. 

As wc see in the following section, the position of 
the minimum is shifted towards higher values of Aq 
for small tan/3. In this case, a change of the sign of 
the bilinear Higgs parameter ^ also has a significant 
impact. 



B. Dependence of the Tree— Level Neutrino Mass 
on the other mSUGRA Parameters 

In App. |X1 we discuss in detail how the neutrino 
mass matrix depends on the other mSUGRA param- 
eter besides ^o- Here we summarize the most impor- 
tant effects and illustrate them in Fig. [51 

In Fig. ini we show two dimensional mSUGRA pa- 
rameter scans of the tree-level neutrino mass. The 
other mSUGRA parameters are those of Point I, 
Sec. IIIBl with A233IGUT = 10~^. One scan param- 
eter is always Aq in order to show how the position of 
the minimum, which was described in the last section, 
changes with the other mSUGRA parameters. 

Fig. 6(a) shows the Aq~Mi/2 plane. We can clearly 



see that the position of the neutrino mass minimum 
is at Aq sa 2M1/2 as was concluded above. This il- 
lustrates that varying M^ /2 has a similar effect on the 
running of /i^ -j,, Eq. ([55)) and Eq. (jSTj), as varying ^q- 
This is clear from the arguments (a)-(c) in Sec. IIV Al 
We could just rephrase the case differentiation as 

(a) Mi/2> Ao/2. 

(b) Ml/2 ~ Ao/2. 

(c) Ml/2 < Ao/2. 
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FIG. 6: Two dimensional plots of the tree-level neutrino mass. The plots are centered around parameter Point I, Sec. Ill Bl 
with A233IGUT = 10"^. The yellow regions signify parts of the parameter space where the neutrino mass becomes smaller 
than lO'^'eV, Figs. [6(a)l [6(c)| [6(3)1 °^ smaller than 10"*eV, Fig.leCb)]. 



For A e {^ijk} the relation is altered to ^0 ~ 
Mi/2/2. The change of the prefactor is due to the 
fact that Xijk couples only leptonic fields to each 
other. Consequently, only superfields carrying SU(2) 
and U(l) charges, but not SU(3) charges, contribute 
to the relevant RGEs; sec Add. lA 5] for more details. 

In Fig. |6(b)] we present the tree-level neutrino mass 
as a function of ^0 and Mq. We observe that the posi- 
tion of the neutrino mass minimum is fairly insensitive 
to Mo, compared to Aq, M1/2 and tan/3 (see below). 
The minimum is shifted to slightly higher values of 
Aq for large Mq. However at large Mq, the interval 
around the minimum in the Aq direction where the 
the tree-level neutrino mass is considerably reduced 
(and therefore the bounds on X'^-j^ are substantially 
weakened) is significantly broadened. 

Finally we show in Fig. 6(c) and Fig. |6(d)"] the Aq- 
tan/3 plane for sgn{ fi) = +1 and sgn(/Lt) = —1, respec- 
tively. In Fig. 6(c) we can see that for low tan/3, the 



This is due to a decrease of the down-type Yukawa 
coupling for low tan/3 leading to a decrease of the RHS 
of Eq. ([55)) . This decrease needs to be balanced by in- 
creasing Aq; recall that h'^j^ = A-^j, x ^0 at Mqut in 
Eq. — 



neutrino mass minimum shifts to higher values of Aq . 



A comparison of Fig. 6(c) and Fig. 6(d) also shows 
that there is a mirror effect around Aq = 800 GcV 
(w 2M1/2) when wc change the sign of the bilinear 
Higgs parameter ^. This happens because of a reversal 
of the sign of the RGE for A, cf. Add. IA3I Therefore 
the shift of the minimum for low tan/3 now appears 
towards lower values of Aq . 



C. The Dependence of the Loop Contributions 

to the Neutrino Mass on the mSUGRA 

Parameters 

The loop contributions to the neutrino mass matrix 
are usually several orders of magnitude smaller than 
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the loop contributions on the mSUGRA parameters. 
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FIG. 7: Ao dependence of the different contributions to the 
neutrino mass at the REWSB scale for the B3 mSUGRA 
Point I, Sec. Ill Bl with A233IGUT = 10"^ Note that only 
the absolute values of the contributions to the neutrino 
mass are displayed. m^J^" and m^'^ are negative whereas 
ml^'' is mostly positive, mj^" is only negative between the 
two minima of \m'f/\: see Sec. IA6lfor details. 
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-loops: 



This contribution to the 
neutrino mass, m'l^'^, depends only weakly on 
the mSUGRA parameters, in particular it de- 
pends logarithmically on the relevant sfermion 
mass. For example, varying Aq from to 1400 
GeV (-200 GeV to 1000 GeV) around Point I 
(Point II) l eave s the magnitude of m^^ nearly 
unchanged |123{ ; cf. the dotted-dashed blue line 
in Fig. [7] (Fig. |8]). However, increasing Mg or 
Ml/2 results in a decreasing m^^: as the SUSY 
spectrum gets heavier the sfermions in the loops 
decouple. 

Neutral scalar-neutralino-loops: This contribu- 
tion to the neutrino mass, m'^'^ , as a function of 
Aq possesses a minimum which lies in the vicin- 
ity of the m}J''° minimum. However, there is no 
exact alignment. This behavior can be under- 
stood by noting that the minima of m^'^ arise 
due to the vanishing of Di, because roughly 



(xDf 



(38) 



cf. App. IA61 This can be seen in Fig. [7]as well as 
in Fig. |5] (dotted magenta line). Again, increas- 
ing Mq or Ml/2 will in general decrease m 
because the SUSY mass spectrum gets heavier. 
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• NLO corrections to the sneutrino vcvs arc typ- 
ically at least one order of magnitude smaller 
than the tree level quantities determining the 
sneutrino vevs, m^ - x VdjiMyYa and Di x 
Vul{Mo)\, in Eq. dfO]) ^^. For illustration, one 
could consider this as a 0(10%) correction to 
m? =: . This shift upwards of the dotted-dashed 
blue line in Fig. |3] slightly changes the position 
of the tree-level neutrino mass minimum, but 
does not alter any of the conclusions drawn in 
this section. Since the effects that we investigate 
in this paper arise mainly from the contribution 
DiVu to the sneutrino vevs (see Sec. lIV A| . these 
corrections are not important for the qualitative 
analysis. 



-200 200 400 600 800 1000 

Aq (GeV) 

FIG. 8: Same as Fig.[7l but for the B3 mSUGRA Point II, 
Sec.HTBl with A233IGUT = lO"'*. 



For parameter Points I and II, Sec. IIIBl the Aq in- 
terval, AAq, where the loops dominate is relatively 
small, c/. Fig. [7] and Fig. [5] However, there are other 
parameter regions where the loops dominate in in- 
tervals of AAq = 0(100 GeV)! This is for example 
the case if one varies Aq around the benchmark point 
SPSla jl04| . We investigate the resulting bounds on 
the LNV trilinear couplings in the following section. 



the tree-level contribution [4l|, 152]. However, in the 
region around the tree-level neutrino mass minimum, 
the loops dominate as shown in Fig. [7] and Fig. [51 
Therefore, we now briefly discuss the dependence of 



V. BOUNDS ON TRILINEAR B3 
COUPLINGS FROM !/-MASSES 

In this section, we calculate upper bounds on all tri- 
linear LNV couplings A e {-^yfc, A'^-^,} at A'/qut from 



13 



the cosiTLological upper bound on the sum of neutrino 
masses as given in Eq. (|4]). Note that in good approx- 
imation 



m„\-EW « A^Igut , 



(39) 



as explained in Scc. l!VAl|T24} . Eq. ^ and App.Ell 
Based on this approximation we employ an iterative 
procedure to account for effects beyond Eq. ((39| . 

In Sec. lv A| we first compare our bounds with those 
given in Ref. [52|, where the mSU GRA parameters of 
the benchmark point SPSla |l04{ (in addition to A) 
were used. We choose the same mSUGRA parameters 
beside Ao in order to show how the bounds change 
in the vicinity of the tree-level neutrino mass mini- 
mum, cf. Sec. IIV Al We then perform in Sec. IV CI 
two dimensional parameter scans around the bench- 
mark scenarios Point I and Point II {cf. Sec. IIIBI) to 
show more generally how the bounds depend on the 
B3 mSUGRA parameters. 

In our parameter scans we exclude parameter re- 
gions where a tachyon occurs [S^l or where the LEP2 
exclusion bound on the light SSM Higgs mass is not 
fulfilled [9II9II. However, we reduce the LEP2 bound 
by 3 GeV in order to accou nt for numerical uncertain- 
tics of SOFTSUSY pMll07| . For instance, in the de- 
coupling limit (where the light Higgs, ft,°, is SM-like) 
a lower bound of 



m,,o > 111.4 GeV 



(40) 



is imposed. In the figures, we also show contour lines 
for the 2a window of the SUSY contribution t o th e 
ano malous magnetic moment of the muon [QJ, llOS 



8.6 X 10- 






10 



(41) 



For more details see Ref. [Gg and references therein. 

We also note that the complete parameter space 
which we investigate in the following (having rejected 
parameter regions which contain tachyons or violate 
the LEP2 Higgs bound) is consistent with the ex- 
perimental upper bound on the branching ratio of 



BR(B,, -^ A^+M") < 4.7 x 10" 



(42) 



and with t he 2a window for the branching ratio of 

6^57, [H [ml. 



2.74 X lO"'' < BR(6 -^ 57) < 4.30 x 10^ 



(43) 



We employed micrOMEGAs2.2 [112| . for the evaluation 
of (5aSUSY^ BR{Bs -^ A'^/""), and BR(6 -> 57). Note 
that there is a significant correlation in mSUGRA 
models between the muon anomalous magnetic mo- 
ment and Bs — > fJ'^l^~ [ll3|. Furthermore, we are well 
above the standard supersymmetric mass bounds, as 
for example on the charginos. 



A. Comparison with Previous Results 

In Ref. [54] , bounds on single couplings A at Mqut 
in the B3 mSUGRA model were determined for the 
mSUGRA parameters of SPSla, in particular Aq = 
— 100 GeV. However, the possibility of obtaining much 
weaker bounds on the coupling A in the region of 
the tree-level neutrino mass minimum was not ex- 
ploited. Note that the bounds in Ref. [53 were also 
obtained for a le ss re strictive cosmological bound of 
^rn,^. < 0.71 cV |l25| . We present here an update of 
these results by using Eq. (|3]). We then explore the 
mSUGRA parameter dependence of the bounds. 

In Tab. Ill and Tab. U (A e {A^^,} and A e {X,jk}, 
respectively), we compare the previous results with 
bounds (at Afcux) that we obtain for identical B3 
mSUGRA parameter points, where only the choice of 
Aq differs. In order to obtain corresponding bounds 
at A/ew one needs to take into account the RGE evo- 
lution of the couplings. Quantitatively this results in 
multiplying the bounds in Tab. | I | (Tab. pB by roughly 
a factor of 3.5 (1.5), cf. Ref. Hi [Hi, Hill le^ . 

In addition to Aq = -100 GeV (SPSla), we choose 
two parameter points which lie AAq w 10 GeV and 
AAq ss 60—70 GeV, away from the neutrino mass 
minimum. In Tab.|T|(A g {A^ j.}), we choose Aq = 500 
GeV (column 3 and 6) and ^0 ~ 550 GeV (column 4 
and 7). In Tab. |lT| (A e {Xijk}), we choose Aq = 200 
GeV (column 3) and Aq = 120 GeV (column 4). This 
enables us to examine the dependence of the bounds 
on Aq around the tree-level mass minimum. 

Note that at SPSla and when varying Aq, the 
neutrino mass minimum for A^ ,, 7^ lies at Aq = 
563 GcV. This value is mostly independent of the 
choice of the indices i,j,k. This is clear because the 
condition for the minimum to occur, Aq w 2Mii2, 
docs not depend on i,j,k, cf. Sec. IIVI Similarly, 
for XijklcvT 7^ the minimum is expected at Aq sa 
Ml/2/2. For t he S PSla parameters we thus obtain 
Ao«127GeV fl26| . 

We first concentrate on Tab. |I| Comparing the 
columns for Aq = —100 GeV and then for Aq = 500 
GeV, i.e. approaching the minimum up to AAq ~ 63 
GeV, the bounds from too large neutrino masses are 
weakened by a factor of 13-15. When we go even 
closer, i.e. Aq = 550 GcV and AAq = 13 GcV, the 
bounds arc weakened by a factor of 40-64 compared 
to Aq = —100 GeV. As we discuss below, in the case of 
up-mixing, some couplings in Tab. U (column 2-4) can 
not be restricted at all by too large neutrino masses. 
In this case we show the bounds at A/qut [marked by 
()*], that one obtains from the absence of tachyons; 
see also Ref. |52| . 

We differentiate in Tab. U between up- and down- 
type quark mixing, cf. Sec. IIIDI Different quark mix- 
ing has important consequences for the bounds on the 
couplings A^ J. if j ^ k. As is clear from Sec. (|IVAp . 
the tree-level neutrino mass is generated proportional 
to A^ -J, X {Y D)jk- Thus, no tree-level mass is gener- 
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Up mixing 


Down mixing 


^0 (GeV) 


-100 


500 


550 


-100 


500 


550 


-^111 


2.0x10"^ 


2.7x10-2 


8.3x10-2 


9.7x10--* 


1.3x10-2 


5.3x10-2 


-^211 


2.0x10-3 


2.7x10-2 


8.3x10-2 


9.7x10-* 


1.4x10-2 


5.3x10-2 


■^311 


2.0xl0--^ 


2.7x10-2 


8.3x10-2 


9.6x10-* 


1.3x10-2 


5.3x10-2 


'^121, Aii2 


(1.3x10-^)* 


(1.7x10-*)' 


(1.7x10-*)* 


4.3x10-* 


6.0x10-3 


2.7x10-2 


'^221, A212 


(1.3x10-^)* 


(1.7x10-*)' 


(1.7x10-*)* 


4.3x10-* 


6.0x10-3 


2.7x10-2 


'^3211 -^312 


(1.3x10^^)' 


(1.7x10-*)' 


(1.7x10-*)* 


4.3x10-* 


5.9x10-3 


2.6x10-2 


•^131 


(1.4x10-^)* 


(1.9x10-*)' 


(1.9x10-*)* 


6.9x10"* 


9.5x10-3 


4.2x10-2 


A23I 


(1.4x10-')' 


(1.9x10-*)' 


(1.9x10-*)* 


6.9x10-* 


9.5x10-3 


4.3x10-2 


-^^331 


(1.4x10^^)' 


(1.9x10-*)' 


(1.9x10-*)* 


6.8x10-* 


9.3x10-3 


4.2x10-2 


■^122 


9.1x10"^ 


1.3x10-3 


5.3x10-3 


8.9x10-^ 


1.2x10-3 


5.2x10-3 


A222 


g.lxlO"'^ 


1.3x10-3 


5.3x10-3 


8.9x10-^ 


1.2x10-3 


5.2x10-3 


-^322 


g.oxio"-"^ 


1.3x10-3 


5.3x10-3 


8.8x10-^ 


1.2x10-3 


5.2x10-3 


•^132 


2.4x10-2 


(1.9x10-*)' 


(1.9x10-*)* 


5.8x10-^ 


8.0x10-* 


3.9x10-3 


-^232 


2.4x10^2 


(1.9x10-*)' 


(1.9x10-*)* 


5.8x10-^ 


8.0x10-* 


3.9x10-3 


A332 


2.4x10-2 


(1.9x10-*)' 


(1.9x10-*)* 


5.8x10-^ 


7.9x10-* 


3.8x10-3 


-^113 


4.2x10-3 


5.5x10-2 


1.9x10-* 


6.3x10-* 


8.7x10-3 


3.8x10-2 


-^213 


4.2x10-3 


5.5x10-2 


1.9x10-* 


6.3x10-* 


8.7x10-3 


3.8x10-2 


-^313 


4.2x10-3 


5.4x10-2 


1.7x10-* 


6.2x10-* 


8.6x10-3 


3.7x10-2 


-^123 


5.9x10-^ 


8.7x10-3 


2.4x10-2 


5.3x10-^ 


7.4x10-* 


3.4x10-3 


-^223 


5.9x10-* 


8.7x10-3 


2.4x10-2 


5.3x10-^ 


7.4x10-* 


3.4x10-3 


■^323 


5.8x10-^ 


8.5x10-3 


2.4x10-2 


5.3x10-^ 


7.2x10-* 


3.4x10-3 


-^133 


2.3x10-'' 


3.2x10-^ 


1.3x10--* 


2.3x10-'^ 


3.2x10-^5 


1.3x10-* 


-^233 


2.3x10-^ 


3.2x10-^ 


1.3x10-^* 


2.3x10-^ 


3.2x10-^ 


1.3x10-* 


-^333 


2.3x10-'^ 


3.1x10-^ 


1.3x10-"* 


2.3x10-'^ 


S.lxlO-'^ 


1.3x10-* 



TABLE I: Upper bounds on the trilinear couplin gs A' ,fc, Eq. ([2]), at Mgut for several values of Aq (second row). The 
other mSUGRA parameters are those of SPSla [104| . We assume up-mixing (down- mixing) in column 2-4 (5-7), cf. 
Sec. Ill D] Bounds arising from the absence of tachyons are in parentheses and marked by a superscript t: ()*. 



Ao (GeV) 


-100 


200 


120 


A211 


1.1x10-* 


2.7x10-* 


(7.1x10-*)* 


A311 


1.1x10-* 


2.7x10-* 


(7.1x10-*)* 


A231 


(5.5x10-*)* 


(6.7x10-*)' 


(7.1x10-*)* 


A122 


4.7x10-* 


1.7x10-3 


4.9x10-3 


A322 


4.7x10-* 


1.7x10-3 


4.9x10-3 


A132 


(5.5x10-*)' 


(6.7x10-*)' 


(7.1x10-*)* 


A123 


(5.1x10-*)* 


(6.3x10-*)' 


(6.7x10-*)* 


A133 


2.7x10-5 


1.0x10-* 


2.8x10-* 


A233 


2.7x10-5 


1.0x10-* 


2.8x10-* 



TABLE II: Upper bounds on the trilinea r co uplings Xtjk, at Mgut for different values of Aq (first row). The other 
mSUGRA parameters are those of SPSla [104| . Bounds arising from the absence of tachyons are marked by ()*. 



ated at this level when we consider j ^ k and up-type 
mixing (which implies a diagonal Yd). But, an addi- 
tional A^j,j, coupling will be generated via RGE running 
at lower scales, cf. Ref. [52[- This coupling will still 
generate a tree-level neutrino mass, which is how ever 
suppressed by the additional one-loop effect (127{ . 

This effect can be seen in Tab. HI if we compare for 
example the upper bounds on A223 and A233 for up- 



and down-type quark mixing. The ratio between these 
bounds is roughly 200 in the case of up-type mixing 
whereas there is only one order of magnitude difference 
for down-type mixing. 



In the latter case, the ratio between the A223 and 
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A233 bounds originates mainly from the ratio 

(Yd)23 _ {^C km) 23 
(Yd)33 (Vc_R-m)33 



(44) 



since the tree-level mass is generated via A223 x (Y£))23 
and A233 X (Y£i)33, respectively. 

To conclude, the bounds from the generation of neu- 
trino masses (at least in the case of down-type mixing) 
are usually the strongest bounds on the couplings A^ ,, 
at Mgut ■ As considered in Ref . [52| , they range from 
O{10-^) to 0(10-6) for the parameter point SPSla 
(column 5 in Tab-IJ). However, there is a large window 
around the tree-level neutrino mass minimum, where 
bounds may be obtained that are between one and 
two orders of magnitude weaker than those in Ref. |52| . 
Around the minimum, the couplings arc only bounded 
from above by ©(IQ-^) to ©(lO"*) (c/. column 7 in 
Tab. [I]). Th us, other low energy bounds become com- 
petitive [4i.[5ll[TTl. 

We now discuss in Tab. |TT] the case of a non- 
vanishing coupling Xijk at A/gut- Contrary to Tab.|Tl 
in the case considered in Tab. |TT]thc quark mixing as- 
sumption does not affect the bounds since Ayfc couples 
only to lepton superfields. Due to the antisymmetry 
Kjk ~ —Xjik there are only 9 independent couplings. 

Wc observe in Tab. HIl that lii ^ j ^ k ^ i there are 
no bounds from too large neutrino masses. The only 
bound we obtain stems from the absence of tachyons. 
This is because we assume a diagonal lepton Yukawa 
matrix Y^; as stated in Sec. IIIDI and therefore, only 
coupli ngs of the form A^fefe can generate a neutrino 
mass [l28l |. 

For these couplings, the bounds at Mqut for Aq — 
— 100 GeV (column 2) range from 1.1 x 10"^ (A211 and 
A311) to 2.7 X 10"^ (A133 and A233). If we approach the 
tree-level mass minimum, i. e. going from column 2 to 
column 4 with Aq = 120 GeV, the bound is weaker 
than the tachyon bound (A211 and A311) or it is weak- 
ened to 2.8 X 10^'' (A133 and A233). The bounds from 
neutrino masses are thus decreased by roughly a factor 
of 10. 

Comparing the bounds on Xikk at Mgut, one can 
see nicely how the choice of k influences the strength 
of the bound. The bounds resemble the hierarchy 
between the lepton Yukawa couplings (Y^;)^^ anal- 
ogously to Eq. (HI]). Therefore, the bounds are 
strongest for fc = 3. 

In contrast to Tab. [H the bounds arc only reduced 
by one order of magnitude when wc approach the tree- 
level mass minimum. This is because the loop con- 
tributions play an important role for the bounds in 
Tab. ini as we discuss in the following section. 



B. Influence of Loop Contributions 

Wc now shortly discuss the influence of the neutrino 
mass loop contributions on the bounds. Typically, one 



expects that the closer we approach the tree-level neu- 
trino mass minimum the more important the loop con- 
tributions become. This is because the loops are not 
aligned to the tree-level mass, cf. Sec. IIVCI 

However, in the case of the neutral scalar loops there 
is still partial alignment, because both the tree-level 
mass minimum and the minima of the neutral scalar 
loops crucially depend on the vanishing of the bilin- 
ear LNV parameter Di, cf. Sect. IIVCI Therefore, 
it is the A'A'-loops and AA-loops, Sec. IIII 21 that are 
relevant whenever the loop contributions become dom- 
inant over the tree-level contributions. 

Wc now give a few examples. For A e {Kjk}, 
Tab.|TTl the loop contributions dominate over the tree- 
level mass in a range of A^o ~ i50 GeV around the 
tree-level mass minimum at ^0 — 127 GeV. Therefore, 
the bounds in this region are much more restrictive 
(i.e. the value of the bounds decreases) when taking 
into account the loop contributions. For example. 



\tot 
^233 

A tree 
233 



0.3. 



(45) 



for Ao = 120 GeV; column 4 in Tab.lHl Here, A*2^*3 is 
the bound on A233 at Mgut if we take into account 
both tree-level and loop-contributions to the neutrino 
mass. In contrast, A23f would be the bound if we only 
employ the tree-level mass. 

Further away from the minimum, the influence of 
the loop contributions is weaker. The bounds are 
strengthened by approximately 5% for Aq = 200 GeV 
(column 3 of Tab. ^ and < 1% for ^0 = -100 GeV 
(column 2 of Tab. Ill- 

The loop contributions arc less important for the 
bounds in Tab. HI i.e. A S {A^-^.}. For example, even 
near the tree-level mass minimum (column 4 and 7 
with Aq = 550 GeV), the bounds become only stronger 
by up to 20% if we take the loop induced neutrino 
masses in addition to the tree-level mass into account. 



Dependence of Bounds on B3 mSUGRA 
Parameters 



In this section, we discuss the dependence of the 
bounds on A G {A^fc, A'-^j,} at Mgut on the B3 
mSUGRA parameters. For that purpose we perform 
two-dimensional parameter scans around the bench- 
mark scenarios. Point I and Point II, of Sec. IIIBI For 
the calculation of the bounds all contributions to the 
neutrino mass considered in Sec. Illll arc included. We 
will focus here on the couplings A233 and A233, be- 
cause these couplings have the strongest constraints 
from neutrino masses, cf. Tab. Hand Tab.lHl 

We have analyzed in Sec. IIVI how the neutrino 
mass changes with the mSUGRA parameters. Due 
to its approximate proportionality to A^, cf. Eq. (jSH), 
the analysis in Sec. IIVI is directly transferable to the 
mSUGRA dependence of bounds on the LNV trilinear 
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(a)Upper bounds on A233 at A/qut in the Aq-Mj /2 plane. 
The parameter space below the green hne is disfavored by 
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(b)Upper bounds on A233 at AfouT in the Aq— tan /3 plane. 

The parameter space below the green line is disfavored by 

Sl^^^; see Eq. glj. 
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(c)Upper bounds on A233 at Mqut in the Aq—Mq plane. 
The parameter space above the green line is disfavored by 
SfVSY. ggg Eq gg 

FIG. 9: Upper bounds on A233 at Mgut from the cosmological bound on the sum of neutrino masses, Eq. Q, as a 
function of mSUGRA parameters. The parameter scans are centered around the benchmark Point I, cf. Sec. IIIBI The 
blackened-out region denotes parameter points where tachyons occur or where the LEP2 Higgs bound is violated. 



couplings. Therefore, the parameter scans presented 
in this section, i.e. Fig. Inland Fig.[TUl resemble closely 
those in Fig. M Sec. El 

We show in Fig. [HI [Fig. [TU] how the bounds on A233 
[A233] at A/gut vary with mSUGRA parameters. We 

[I0(a)}[l0(c)] the Ao 



present in Figs. 9(a) 9(c) 



M 



1/2 



Ao-tan/3, and Xq 



Figs. 
A/o planes, respectively. The 



bounds are shown on a logarithmic scale. The black- 
ened out regions designate areas of parameter space 
which are rejected due to tachyons in the model or vi- 
olation of the LEP2 bound on the lightest Higgs mass, 
cf. Eq. (Uni)- Furthermore, we include contour lines 
of the 2(T window for the SUSY contribution to the 
anomalous magnetic moment of the muon, Eq. (|4ip. 
Imposing Eq. (pij) disfavors the parameter space be- 
low [a bove the g reen contour line i n Fig s. 9(a)[ |9(b)[ 
10(a)| and [TO(b)] [Fig. f9(c)] and Fig. [9(c)] . 

We observe in Fig. [9] that the strictest bounds on 



A233 from too large neutrino masses are of 0(10^^) . 
However, there are sizable regions of parameter space 
where the bounds are considerably weakened. For ex- 
ample, in the A0-M1/2 plane. Fig. 9(a)[ the bounds 
are of 0(10"^) only in approximately half of the pa- 
rameter space whereas in the other half, the bounds 
are 0(10"^) or weaker. In roughly 10% of the al- 
lowed region in Fig. [SJ the bounds even lie at or above 
0(10"^)! In this region, the loop contributions to the 
heaviest neutrino mass are essential for determining 
the bounds since the corresponding tree-level neutrino 
mass vanishes, cf. also the discussion in Sec. IV Bl 

We can see in Fig. [10] a similar behavior for the pa- 
rameter dependence of the bounds on A233. Here, the 
strongest bounds arc now of 0(10^''). However, for 



example in the Aq-Mq plane. Fig. 10(c) the bounds 
are as strong as 0(10^^) in only about 25% of the 
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FIG. 10: Same as Fig. |9l but for A233 at AfcuT and for the benchmark scenario Point II, cf. Sec. IIIBI 



parameter plane. The remaining 75% have bounds of 
0(10-*) (50%) or even 0(10-3) (25%)! 

Up to now, we have analyzed how the bounds on 
the trilinear LNV couplings A233 and A233 vary with 
the mSUGRA parameters. However, from the analysis 
in Sec. IV Al we can easily deduce how most of these 
bounds change for different couplings A'j ■;, and A, 



\jk 1 



for different indices i,j,k. For A^ j, the index i does 
not significantly influence the bound, because the em- 
ployed Yukawa coupling, (Y£))jfc, via which the tree- 
level mass is generated, does not depend on i. But, the 
situation is totally different when wc change the in- 
dices j, k. In general, for A' ;, (and down-mixing) the 
bounds will display the hierarchy of the down-type 
Yukawa couplings. Therefore, bounds for couplings 
X^■^^■^ arc about three orders of magnitude weaker than 
bounds for the couplings A^33 as long as the other B3 
mSUGRA parameter are the same. We also observe a 
similar behavior for A^ j, with up-mixing and for Xijk 
[using {YE)jk instead of {YD)jk], if j = k; cf. the 
discussion in Sec. IV Al 

To conclude, one can use the Yukawa matrix Yjj 
(Ye) to easily translate the bounds in Fig.l9l(Fig. [T0)) 



to bounds on couplings other than A233 (A233) 



VI. SUMMARY AND CONCLUSION 

We have calculated upper bounds on all trilinear 
lepton number violating couplings at the grand unifi- 
cation scale within the B3 (i.e. lepton number violat- 
ing) minimal supergravity (mSUGRA) model, which 
result from the cosmological bound on the sum of neu- 
trino masses. We have shown that these bounds on the 
couplings can be weaker by one to two orders of mag- 
nitude compared to the ones which were previously 
presented in the literature for the benchmark scenario 
SPSla; cf. Sec. |Vl In general, the bounds can be 
as weak as 0(10"^). Thus other low energy bounds 
become competitive. 

The reason for these large effects is that the tree- 
level neutrino mass depends strongly on the trilin- 
ear soft-breaking Ao-parameter (and also similarly 
on the gaugino masses). We concluded in Sec. ITVl 
that in regions of parameter space with Aq w 2Mi/2 
[Ao « A/1/2/2) for A^Jgut 7^ (A,,fc|GUT 7^ 0), a 
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cancellation between the different contributions to the 
tree-level mass can occur. We have explained this ef- 
fect in detail and have shown that such a cancellation 
is significant in large regions of the mSUGRA param- 
eter space. For example, the bounds can be weakened 
by one order of magnitude in Ao intervals of up to 
0(100 GeV), see Figs. M and [TOl Therefore, much 
weaker bounds (compared to previous ones) can be 
obtained without significant fine-tuning. 

In order to obtain the correct bounds in the vicinity 
of the tree-level neutrino mass minimum, we included 
the main loop contributions to the neutrino mass ma- 
trix; cf. Sec. mil We also described in Sec. IIVI and 
App. |A]for the first time the dependence of the tree- 
level and loop induced neutrino mass on all mSUGRA 
parameters. Although we concentrated in this work 
on the B3 mSUGRA model, the mechanisms described 
will also work in more general R-parity violating mod- 
els. 

Our work can help to find new supcrsymmetric sce- 
narios that are consistent with the observed neutrino 
masses and mixings. We have shown in this publica- 
tion how the (typically large) hierarchy between the 
tree-level and 1-loop neutrino masses can systemat- 
ically be reduced. Together with at least one addi- 
tional lepton number violating coupling, one can use 
this mechanism to match the ratio between tree-level 
and 1-loop induced masses to the observed neutrino 
mass hierarchy, both for hierarchical neutrino masses 
and for a degenerate spectrum. 

Wc also note, as described in the introduction, that 
large lepton number violating couplings can lead to 
distinct collider signatures. We will address these top- 
ics in future publications. 
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Appendix A: ;/-Masses: Dependence on Further 
B3 mSUGRA Parameters 



In Sec. IIVAI we described in detail the dependence 
of the tree-level neutrino mass, Eq. ([T7)) . on the B3 
mSUGRA parameter ^o- We also reviewed some fur- 
ther effects in Sec. IIVBI In this appendix, we explain 
now in more detail the dependence of the tree-level 
neutrino mass and the loop induced masses on the re- 
maining B3 mSUGRA parameters. 




100 200 300 400 500 600 700 800 900 1000 

Ml/2 (GeV) 

FIG. 11: Same as Fig. [1 but now for the mSUGRA pa- 
rameter Ml/2 instead of Aq. 



1. Afi/2 Dependence 

The tree-level neutrino mass minimum can be ex- 
plained equivalently in terms of its dependence on 
Afi/2 instead of its dependence on ^o- This is because 
varying M1/2 has a similar effect on the sneutrino vev 
Ik, Eq. dSJ), as varying Aq, cf. Sec.|lVA]and Scc. lTVBl 
However, when varying M1/2 there are additional ef- 
fects coming on the one hand from the dependence 
of /i^, {M?)ii and m^ - on Mi/2- These quantities 

are linear functions of My^ . For //^ this can bee seen 
from Eq. ([M|) . For {M?)u and m^ - this follows be- 
cause the respective RGEs are functions of the squared 
sfermion masses [52| . One obtains for example 



(M|). 



M^ + 0.52M 



1/2 



1 9 

-Mj cos 2/3. 



(Al) 



On the other hand, there is also a direct proportional- 
ity of m|^'™ to M~\, cf Eq. HH). All these additional 
effects do not significantly influence the position of the 
tree-level neutrino mass minimum, i.e. Aq « 2Mi/2 
still holds for A e {A^-^.}; see Sec. |lVl However, the 
effects add a global slope to the terms (as a function of 
-^^1/2)1 which contribute to the tree level mass. This 
behavior can be seen in Fig. 1111 

Wc show in Fig. [TT] the same contributions as in 
Fig. [31 but now as a function of M1/2 instead of Aq. 
Here Aq has been fixed to 900 GeV. On the one hand, 
we observe that the quantities DiVu (dotted magenta 
line) and (m^ - + fiKiJvd (dottcd-dashcd blue line) 
are nearly constant for low values of M1/2, but they 
have a positive slope for large values of Mim. This 
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FIG. 12: Same as Fig. E] but now for the mSUGRA pa- 
rameter tan j3 instead of Aq . 



is mainly due to their dependence on /i; c/. Eq. ^ 
[Eq. (g])] for bi [ni\. On the other hand {M^'^)ii (sohd 
turquoise hnc) has a negative slope for all values of 
Ml/2 because of Eq. (jAip . Overall this leads to a steep 
decrease of the tree-level neutrino mass (solid red line) 
in the region of low Mi/2, whereas in the region of large 
Ml/2, the various contributions' dependence on M1/2 
roughly cancels, see Fig. [Til 

Going beyond the plot, for M1/2 — > 00 the tree-level 



mass scales with M. 



1/2- 
tree 



as follows from the different 



contributions to vt.^J'"^ in Eq. (|17p . Such a behavior is 
expected, because SUSY decouples from the SM sector 
in the limit Mi 12 — > cxd. 



2. tan/3 Dependence 



This (weak) tan/3 dependence of ll^iWul is illustrated 
in Fig. [12] for our B3 mSUGRA parameter set Point I; 
see Sec. IIIBI One observes that the dotted magenta 
line {\DiVu\) increases between tan/3 = 2 and tan/3 w 
40. Here, sgn(c2) > 0. Above tan/3 « 40, |Aw„| 
starts decreasing, i.e. sgn(c2) < 0. This is due to 
the enhancement of the down-type Yukawa coupling 
when increasing tan/3, since this reduces /i^j, further 
and further until it becomes negative. This decrease 
of |i^iZ^„| is only partially visible in Fig. [T^ since the 
parameter region with high tan^ is excluded due to 
tachyons. 

One can also see in Fig. [12] that the other term de- 
termining the sneutrino vev, (m^ - -I- ^kAv^,, which 
is displayed as a dotted-dashed blue line, is fairly con- 
stant regarding tan/3. This contribution to the sneu- 
trino vev is subtracted from the first term, DiVu (dot- 
ted magenta line), so that the sneutrino vev becomes 
zero when the two lines intersect; see Eq. (|?T|) . 

We observe this intersection in Fig.[T2]at tan/3 « 22, 
thus yielding the tree-level neutrino mass minimum 
in this region. In theory, there could even arise two 
minima because above tan/3 sa 40 DiVu starts de- 
creasing again, leading to another intersection with 
(m^ - -|- ^Ki)vd. However, as mentioned before, this 
usually happens in an excluded region of parameter 
space. 

As is also illustrated in Fig. [T^J there is quite a siz- 
able difference between the two terms which determine 
the sneutrino vev, ie. (mj = +/iK,;)ud (dotted-dashed 

blue line) and DiVu (dotted magenta line) in the re- 
gion of low tan /3. If we are looking for a neutrino mass 
minimum in this region of parameter space, we need 
to adjust Aq towards higher values, which will increase 
^'ijk I'^f- ^1- (EHj)]- Therefore, increasing Aq will shift 
the dotted magenta line upwards until it intersects 
with the dotted-dashed blue line at the desired low 
tan/3 value. This shift of the tree-level neutrino mass 



minimum to higher Aq is clearly visible in Fig. 6(c) 
For tan/3 = 20, the minimum lies at Aq « 900 GeV 
whereas for tan/3 = 5, it has shifted to Aq « 1300 
GeV. 



Varying tan/3 most importantly affects the tree- 
level neutrino mass via the term DiVu in Eq. pip . The 
RGE for Di, Eq. (|35]) . is proportional to the down- 
type Yukawa coupling {YD)jk = {'m.djjklvd- There- 
fore, 



DiVu ex Ci +C2 — 

Vd 



ci + C2 tan /3 



(A2) 



at A/ew. The factors ci and C2 depend on the other 
mSUGRA parameters but their magnitude is approx- 
imately independent of tan/3. However, there is a de- 
pendence of sgn(c2) on tan/3 via the RGE of h[-^. Es- 
pecially in case (b) of Sec. IIV Al i.e. in the region 
around the tree-l evel neutrino mass minimum, this 
becomes relevant ll29l. 



3. sgn(^) Dependence 

A change of sgn(/i) notably affects the tree-level 
neutrino mass via the RGE running of Di [Eq. (P5|) ] , in 
which the overall sign of the RGE is changed. There- 
fore, the sign of Di itself is reversed at any energy scale 
but its magnitude is mostly unaffected. Consequently, 
the Aq value where 1?^ = is still mostly the same 
after a sign change. 

However, at the position of the tree-level neutrino 
mass minimum, Di needs to be slightly larger than 
zero in order to cancel the other terms contributing to 
the tree-level mass, c/. Sec. IIV Al and Sec. lA2l When 
we are at a parameter point where the tree-level neu- 
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trino mass mininium occurs for positive fj. (i.e. Di is 
small and positive) , a sign change to sgn(/i) = — 1 will 
yield a Di which is small and negative. The other con- 
tributing terms undergo no overall sign change. If we 
would like to obtain a neutrino mass minimum now, 
Di needs to be increased in order to become slightly 
larger than zero again. This can be achieved by de- 
creasing Aq, Sec. IIV Al (or, equivalently, increasing 
Ml/2, Sec. lA 1[) since this increases Di via ft., ;, in its 
RGE, Eq. psp . when /i is negative. Therefore, the 
tree-level minimum will occur at smaller values of Aq 
(or equivalently larger values of M1/2) when we change 
sgn(/^) = +1 to sgn(/i) = -1. 

This effect becomes more important when we go to 
regions of low tan/3. Here the influence of /i^ j. on 

Di, Eq. ((35)) . becomes weaker due to the decrease of 
the down-type Yukawa coupling, as we discussed in 
Sec. IA2I In order to still obtain a positive Di after 
reversing sgn(//), h'^-^. has to decrease in a more sub- 
stantial fashion than for large tan f3. Therefore, the pa- 
rameter point where the tree-level neutrino mass min- 
imum is located will shift to smaller Aq when changing 
sgn(/x) = -|-1 to sgn(^) = —1, especially for tan/3 < 10. 

Overall, this leads to a "mirroring" of the tree-level 
mass minimum curve in the Ao~tan/3 plane around 
Ao = 80 GeV(« 2M1/2). This can be seen in 
Fig. 6(c) and Fig. |6(d)[ for sgn(/i) = +1 the min- 
imum shifts to higher values of ^0 with decreasing 
tan/3, whereas for sgn(/i) = — 1 the minimum shifts to 
lower values of Aq. 



4. Mo Dependence 

Varying Afg does not greatly affect the tree-level 
neutrino mass. However, similar effects as those de- 
scribed in Sec. lA II as additional effects, arise due to 
the dependence of several parameters on Mq , cf. for 
example Eq. (|M)) and Eq. (JA1[) . This can be seen in 
Fig. 1131 where we again show the terms, which enter 
the tree-level neutrino mass formula, Eq. pT)) . We can 
see that most of the quantities depend only weakly on 
Mq. This results in a nearly constant tree-level neu- 
trino mass, cf. solid red line in Fig. [T3l 

However, the above mentioned Mq dependences 
lead to a moderate shift of the tree-level neutrino mass 
minimum towards higher values of ^0 when increasing 
Mq. Explaining this in detail is fairly lengthy because 
the Mq dependence of the parameters determining the 
tree-level neutrino mass is not as straightforward as 
the dependence on other mSUGRA parameters. How- 
ever, the effect is shown numerically in Fig. |6(b)[ 

It should be noted that there is a similar mirror 
effect when changing sgn{fi) as for tan/3. For sgn{fi) ~ 
— 1, the minimum shifts towards lower values of ^0 
when increasing Mq. 
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FIG. 13: Same as Fig. H but now for the mSUGRA pa- 
rameter Mo instead of Aq . 
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FIG. 14: Sa me as Fig. [3 but now for the B3 mSUGRA 
Point II, Sec UlB] with A233IGUT = lO""*. 



5. Changes for A G A 



ijk 



We now consider the case of A G {A^fe} instead 
of A e {Kjk}- Since \ijk only couples lepton super- 
fields to each other (as opposed to the A'- j. operator 
which also involves quark superfields) , the RGEs in 
Sec. IIV Al arc reduced by a (color) factor of 3 [5al54|. 
In addition, the down quark Yukawa matrix elements, 
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(Y]j)jk, need to be replaced by the respective lep- 
ton Yukawa matrix elements, {YE)jk- Otherwise, the 
structure of the RGEs remains the same. 

The only RGE where there arc more extensive rele- 
vant changes is that for hijk (which replaces /i^.j,); cf. 
Eq. dS]). Eq. dSH) must be replaced by [H 



167r^ 



,dh 



ijk 



dt 



9 

-gf{2MiXijk - Kjk) 



+3<7|(2M2Ayfc - h. 



ijk) 



-(A3) 



with hijk = ^0 X Kjk at A'/qut- This looks exactly 
the same as the RGE for h'--/^, Eq. ([55)) . only with 53 
and M3 replaced by ga and Ma {a = 1,2). However, 
it is important to realize that the running of ga and 
Ma is different from the running of 53 and M^. As 
was mentioned in Sec. IIV Al the latter quantities in- 
crease when runnin g to lower energy scales whereas 
the former decrease [23| . 

This has important consequences for the position 
of the tree-level neutrino mass minimum. The terms 
gaMaXijk of Eq. (jA3[) now decrease [as opposed to 
g^M^X'-jf. in Eq. ([55]) ]. It is thus necessary to choose 
Aq smaller in order to have a smaller hijk at AfcuT 
and at lower scales to compensate for this. Quan- 
titatively, we checked numerically that we now need 
Aq « Ml/2/2 (A S {Kjk}) to achieve a vanishing 
tree-level neutrino mass rather than Aq w 2Mii2 
(A g {X'^jk}) as was the case in Sec. IIV Al 

For illustrative purpose, we show in Fig. [14] the Aq 
dependence of the tree-level neutrino mass (solid red 
line) and of the terms determining the sneutrino vev V2 
for a non- vanishing coupling A233 at Mgut- Fig. [HI is 
equivalent to Fig.[3]bcsidc the fact that we now employ 
the parameter Point II with A233IGUT — 10"^ instead 
of the parameter Point I with A233IGUT — 10"^, cf. 
Sec. IIIBI The qualitative behavior of all terms is the 
same in both figures. However, in Fig. [TJ] the minima 
are shifted to lower values of Aq compared to Fig. [3] 

We conclude that the line of argument explain- 
ing the minimum of the tree-level neutrino mass in 
the case of A e {Kjk} ^^^^^ holds for A G {A^fc}. 
However, the position of the minimum now shifts to 



Aq « A/1/2/2. 



6. Ao Dependence of the Neutral 
Scalar— Neutralino— Loops 

According to Eqs. ([26]) and ([27|) . the dominant loop 
contribution from neutral scalar-neutralino-loops to 
the neutrino mass matrix, {m'/,'^)ii, is proportional to 

(m^'^)ij oc {DiVd - Bv^f 

x/(m|o,m|^,m|^^^,m^^,m^J, (A4) 

where / is a function of the neutralino, sneutrino and 
Higgs masses squared, respectively. 

The Aq dependence of Eq. (|A4p is mainly deter- 
mined by Di , since the Aq dependence of Vi is governed 

by A(Ao), 



u,(Ao)cx A(Ao)-t-c 



(A5) 



where the term c depends mainly on the other 
mSUGRA parameters but barely on Aq, as discussed 
in Sec. IIVAI Therefore (to^'^)^^ is roughly propor- 
tional to Df. The behavior of Di has been discussed 
in detail in Sec. IIVAI in the context of the tree-level 
neutrino mass. We have shown that there is always 
a value of Aq where Di becomes zero. Thus the neu- 
tral scalar-ncutralino loops display a similar minimum 
as the tree-level neutrino mass. The position of the 
minimum is close to the tree-level one, but not ex- 
actly aligned. This can be seen by comparing the dot- 
ted magenta line and dashed green line in Fig. [7] and 
Fig. [5] However, since Eq. (|A4p is only an approxi- 
mate formula [for the exact formula, cf. Eq. ([^^ ]. the 
real curve is slightly shifted downwards such that its 
minimum reaches negative values. Therefore \{fn'/,'^)ii\ 
in Fig. [7] and Fig. [5] appears to have two minima. 

It is also immediately obvious from Eq. (|A4[) that 
the scalar-neutralino-loops are roughly proportional 
to [A X (Y£))jfe]^ like the tree-level mass. 
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From m^J^'^ oc Vi we then obtain the relationship in 

Ref. [52] did not calculate the dominant loop contri- 
butions to the neutrino mass matrix. However, these 
are negligible for SPSla, because SPSla lies far away 
from the tree-level mass minimum. 
This value is smaller than would be expected by 
estimating Aq « Mi/2/2 — 250 GeV, because we 
are considering a parameter point with relatively 
low tan/3 (tan/3 — 10 for SPSla). As discussed in 
Sec. IIVBI this leads to a shift of the tree-level neu- 
trino mass minimum towards lower values of Aq, cf. 



Fig. 6(c 



Note that also the loop contributions are strongly 
suppressed, because the A'A'-loops are proportional 



to A' 



X'ikj, Eq. H18p . and the neutral scalar loops 



are aligned with the tree-level mass, cf. Sec. IA6I 
This would change drastically if the Y_b were 
strongly mixed [52 ]. 

In case (a), C2 remains always negative and in 
case (c), C2 is positive. 



